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A hybrid system composed of an isotropic nanoparticle and a semi- 
conductor heterostructure with a quantum well has been consid- 
ered. The nanoparticle is supposed to be polarizable in an external 
electric field. A theoretical model of the hybrid system is substan- 
tiated and formulated. Exact solutions of the model equations 
are obtained. The frequencies of charge oscillations in the hybrid 
system and their damping owing to the dipole— plasmon interac- 
tion are found, the damping mechanism being similar to that of 
Landau damping. The space-time behavior of concentration per- 
turbations in the two-dimensional electron gas is analyzed, and 
the polarization oscillations of a nanoparticle are studied. The 
induced polarization of a nanoparticle at nonzero electron drift 
velocities is found to have a complicated dynamics. In particu- 
lar, the polarization vector circulates along elliptic trajectories for 
two of three frequency dispersion branches. If the electric current 
flows through the quantum well due to an applied electric field, 
the damping of oscillations in the hybrid system is replaced by 
their growth in time, which corresponds to the electric instability 
of the system. New phenomena in hybrid systems can be used to 
excite the emission of nanoparticles by an electric current and to 
electrically stimulate the emission in the terahertz spectral range. 



1. Introduction 

Physical phenomena in the terahertz (THz) spectral 
range are intensively studied, which is associated with 
the fundamental importance of new terahertz physics, as 
well as with a considerable number of feasible applica- 
tions. In particular, many researches aim at constructing 
new sources and detectors of THz radiation. 

In this context, of great interest are the researches 
of semiconductor heterostructures with quantum wells 
(QWs), in which collective oscillations of a two- 
dimensional electron gas (2DEG), i.e. plasmons, can be 



excited. The frequencies of plasmon oscillations belong 
to the THz spectral range. In a homogeneous 2DEG, 
plasmons are stable, i.e. their oscillations attenuate in 
time owing to the electron scattering by various crystal 
defects and by means of the known mechanism of Lan- 
dau damping. Various methods were proposed in order 
to achieve the instability and obtain oscillation-growth 
effects intended to be used for the amplification and the 
generation of THz radiation. Among those, a possibility 
of the instability excitation owing to the electron drift 
in an electric field [H-0> various variants of two-beam 
instability and others were examined. However, 

for such instabilities to be realized in a spatially uni- 
form plasma, very high drift velocities, which are diffi- 
cult to be obtained experimentally, are required. The 
situation essentially changes in spatially non-uniform or 
finite systems For example, as was shown in works 
[6], if special types of contacts are applied to a quan- 
tum well with a finite length, the electron system be- 
comes unstable already at moderate drift velocities of 
electrons. This instability has been studied experimen- 
tally [7rl|. 

Another class of objects, which are active in the THz 
spectral range, includes quantum dot s fl0rll3j . molecules 
and some molecular compounds shallow impu- 

rity centers 17 - HI , and so forth. For brevity, let us call 
such "zero-dimensional" objects as nanoparticles (NPs) . 

Hybrid systems consisting of nanoparticles and het- 
erostructures with free electrons constitute a new type 
of heterodimensional objects, which should demonstrate 
essentially new properties and effects. In particular, if 
the oscillation frequencies of zero-dimensional NPs and 
two-dimensional plasmons are in the THz spectral range, 



ISSN 2071-0186. Ukr. J. Phys. 2012. Vol. 57, No. 3 



367 



V.A. KOCHELAP , S.M. KUKHTARUK 







Z 



Y 



NP 

h 



QW 



Fig. 1. Schematic drawing of the hybrid system. The OY axis is 
perpendicular to the figure plane 



hybrid systems may reveal new properties at extremely 
high frequencies. Therefore, the study of the interaction 
between such NPs and plasmons under equilibrium and 
nonequilibrium conditions seems well-timed. 

In work [20], the studies of the interaction between 
a remote nanoparticle and a drifting 2DEG at rela- 
tively low frequencies, i.e. when the application of the 
"low-frequency" drift-diffusion approximation is valid, 
were started. In this work, we present the results 
of our researches of THz-frequency properties of hy- 
brid systems in the "high-frequency" approximation for 
the electron motion in a 2DEG. We revealed collec- 
tive oscillations of interacting nanoparticles and the 
2DEG. In particular, we determined the frequencies of 
collective oscillations and their extra damping stem- 
ming from the interaction between the nanoparticle 
and plasmons. The damping character is similar to 
that of Landau damping. This system was demon- 
strated to become unstable under definite nonequilib- 
rium conditions, i.e. there emerges an increase of col- 
lective oscillations in time induced by an electric cur- 
rent. The nanoparticle destroys the translational sym- 
metry, and the system becomes non-uniform, which 
promotes the appearance of an instability, similarly to 
what takes place in non-uniform or spatially confined 
plasma [5|. 

The structure of the paper is as follows. In Section 2, 
the model of a hybrid system is discussed, and the ba- 
sic equations are substantiated. In Section 3, the exact 
solutions of the equations are obtained, and their gen- 
eral properties are examined. In Section 4, the main 
attention is given to the stability/instability of a hy- 
brid system under the condition of two-dimensional elec- 
tron drift, the numerical results are presented, and the 
specific examples of nanoparticles and heterostructures 
with two-dimensional electrons are considered. The fea- 
tures of coordinated space-time charge waves are dealt 
in Section 5. The properties of the NP dipole moment, 
plasmon field, and emission by the hybrid system are 
studied in Section 6. Conclusions are presented in Sec- 
tion 7. 



2. Model and Basic Equations 

Consider a hybrid system represented schematically in 
Fig. 1. It consists of a spherical nanoparticle, which can 
be polarized in an external field, and the heterostructure 
with a quantum well, in which two-dimensional electrons 
are localized. The distance between the nanoparticle and 
the quantum well is h. The electrons and the nanoparti- 
cle interact by means of the electrostatic field E = — V0, 
where cf>(x, y, z, t) is a self-consistent electrostatic poten- 
tial. The latter is described by the Poisson equation 



4tt 



A<f> = e(n - n )6(z) + ( D(i)V) S(x)S(y)S(z - h), (1) 



where k is the dielectric permittivity, e the electron 
charge, n(x, y, t) the surface concentration of electrons, 
no the equilibrium surface concentration of electrons, 
D(i) the electric dipole moment of a nanoparticle, and 
8{x) the Dirac delta- function. 

It is convenient to express the scalar potential as a 
sum of two terms associated with the electrons and the 
NP, (/>(x,y,z,t) = <p e (x,y,z,t) + (j> d (x,y,z,t), where 4> e 
and (pd are determined from the system of equations 



KA<f) e — Ane{n — no)<5(z), 
nA(j) d = 4tt (D(i)V) S(x)S(y)S(z 



h). 



(2) 



To describe electrons in the quantum well (the plane 
z = 0), we use Euler's equation of motion and the con- 
tinuity equation, which couple the surface concentration 
n(x, y, t) and the electron velocity t>(x, y, t) with the elec- 
trostatic field. In the general case, we assume that elec- 
trons move with an average drift velocity vo. Then, the 
equations are Q 



— + (t)V||)0= -V0 
ot m 



- v 



On 



+ V|i (tm) = 0, 



(3) 



where m* is the effective mass of an electron, the term 
°~ Vp makes allowance for the electron scattering by crys- 
tal defects, and t p is the time of electron momentum 
relaxation. The subscript |j is used for quantities and 
differential operators dealing with the plane z = only. 
For instance, the Laplace operator in the plane z = is 
A || = ^3 + ^p-. At the same time, the Laplace operator 

in the volume is A = Ay + ^-y. We will also mark the 
real and imaginary parts of complex numbers as primed 
quantities, e.g., Z = Z' + iZ". 

In the two-dimensional electron gas, there emerge col- 
lective oscillations of the charge density, which are called 
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plasmons. The hydrodynamic model with Eqs. (0) is 
valid, when the frequency of these oscillations, w, and 
the absolute value of the wave vector, k, satisfy the con- 
ditions 



OJTp > 1, ql p ~ lp/h » 1 



(4) 



where l p is the mean free path of electrons. The for- 
mulation of the problem includes only one characteristic 
length scale; it is the distance h between the plane of the 
2DEG and the NP. Therefore, one may expect that the 
wave vectors of plasmons, which make the major con- 
tribution to the interaction between the electrons and 
the NP, should be of the order of 1/h. Note that in- 
equalities ((3]) correspond to the ballistic character of the 
electron motion in the spatial region, which is actual for 
the nanoparticle-plasmon interaction. 

Suppose that the velocity of electrons is 
v(x, y, t) = vo + v(.t, y)e~ lut , and their concentra- 
tion n(x, y, t) = no + n(x, y)e~ lut , with v <C vq and 
n <C no- Let the time dependences of potentials also 
vary harmonically, i.e. <fi e (x,y, z,t) — (p e (x,y, z)e~ lu)t 
and <j>d(x,y, z,t) = tpd(x,y, z)e~ tuJt . The frequency w 
in these formulas should be determined by solving the 
equations; generally speaking, it is a complex- valued 
quantity, which corresponds to the oscillations (damped 
or growing in time) in the hybrid system. (If w" > 0, 
the system is unstable, and uj" is called the instability 
increment.) The quantities \{x,y) and n(x,y) de- 
scribe the spatial dependences of electron velocity and 
concentration perturbations, respectively. 

According to Fig. 1, the nanoparticle is located at the 
point (x = 0, y = 0, z = h). In an external electric field, 
it becomes polarized and can be characterized by an elec- 
tric dipole moment, which is determined by the relation 

D(t) = d(w)e-*"", 



d(w) = - Po(u)VVe\ x=y = 0;Z = h , 



(5) 



where /3o(w) is the nanoparticle polarizability. Note that 
formula ([5]) includes only the electron potential and, 
hence, the dipole self-action is excluded from consider- 
ation. If the nanoparticle is isotropic, its polarizability 
can be expressed in the standard form (2l| 



2uj f e* 



K/Mo)| s 



— UJ. -f ibj/Tf ' 



(6) 



where h is the reduced Planck constant, fiw/ is the en- 
ergy of the /-th nanoparticle level, r/ is the lifetime of 



this level, (/|x|0) is the matrix element of the dipole 
transition between the ground and /-th excited states. 

The largest contribution to polarizability ^ ) is given 
by the transition between the ground and the first ex- 
cited state. In this case, the polarization is 



A)M = — j- 



l(iHQ)l 2 

h uj — uj + 270 



(7) 



where the characteristic frequency of dipole oscillations, 
wo, and the damping of dipole oscillations corresponding 
to the transition between the ground and the first excited 
state, 70 = 1/(2^), are introduced. While deriving for- 
mula ([7]), an assumption was made that ui is close to wo, 
because the main effects are expected to take place at 
frequencies close to the characteristic frequency of dipole 
oscillations. 

Equations ([TJ) — , ([5|), and ([7|) compose a basic system 
for the problem under consideration. 

3. Solutions of the Equations and Their General 
Properties 

The dipole potential, being the solution of the second 
Poisson equation in system ([2]) , has the well-known form 



1 d x (ui)x + d v (u)y + d z (uj)(z — h) 
<pd = - - 

K 



(8) 

(x 2 +y 2 + (z-h) 2 )i 

where d x , d y , and d z are the components of the dipole 
moment d(w). Below, for all quantities depending on 
coordinates, the two-dimensional Fourier transformation 
is used: 



<p e (T,z) = / d 2 kip k (z) e 



n(r) = / d\ n k e ikr , V(r) = / d\ 4> k e 



zkr 



(9) 



(10) 



where r= (x,y), k= (k x ,k y ), and V(r) = ip d (r, z)\ z=Q . 
The quantities ifk(z), nk, and ipk are the Fourier trans- 
forms of the electron potential, surface concentration, 
and potential induced by the dipole in the quantum-well 
plane (z = 0), respectively. Then, the Poisson equation 
for ifk(z) and the required boundary conditions are 



d 2 ^ 
dz 2 

d<p+ 



- = o, 



>±oo 



dz 



^0, 



(11) 



z=0+e 



0-e ' 



z=0- 



Airenk 
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where ip^ (z) 



z < 0, k 



^>k(z) at z > and ip k (z) 
+ k 2 , and e +0. 



ifk(z) at where the notations 



The solution of system (fTTj) and the functions ipk and 

rife are 



27ren 



ipk = 
rik = — 



nk 
1 



A X k X + dyky ^ \ ^_ kj/ 

2ttk V k 



(12) 



enok 2 



m (u) - v k)(u; — v k + i/r p ) 



Note that, in the limiting case where the dipole is in- 
finitely far from the 2DEG (i.e. h — > oo), solutions (|12p 
easily bring about the dispersion law for collective ex- 
citations in the drifting 2DEG, i.e. drifting plasmons, 



uj± = vok ± 



27re 2 



no_ k l p tip 

Km 4 2 



(13) 



where the notation j p = 1 /t p is used. On the other hand, 
the zero of the denominator of /3o(w) (Eq. ([7|)) corre- 
sponds to the frequency and the damping of a dipole, 
lo = oj - «7o- 

If the distance h is finite, system (fT2]) and relations 
<(5]) and dl|) can be used to derive an integral equation, 
e.g., for the Fourier transform of the concentration of 
electrons interacting with the dipole, 



e n (3 (uj)ke ' l (k x I x + k y I y + kl z ) 
K 2 m ( w - Vo k)( W - v k + i lp ) - 



(14) 



where I x , I yi and I z are the functionals of rif., 
I x =Jd 2 q^n q e~o h , I y ^Jd 2 q%i q e^ , \I z ^Jd 2 qn q e 



■qh 



With the notations 
B = —/3 (uj), 



A e (w, k) = (oj - v k)(o; - v k + ij p ) 



2ire 2 nok 



the integral equation (fl"4"|) yields the following system of 
algebraic equations: 



Ix — B (S X I X + SoI z ), 

Iy = Bq Sy Iy , 

Iz = Bq(SqI x + S Z I Z ), 



(15) 



S x = 



s z = 



-2kh 



A e (u,k) 



d 2 kk 2 e~ 2kh 

A e (w,fc) ' 



■ ,S y — 



-2kh 



A e (u>,k) 



S = 



A e (w,fc) 



(16) 



are used for the calculable integrals. Note that the pa- 
rameter So = 0, if the electron drift is absent (vo = 0). 

The zero value of the determinant of the system of 
equations (|15p is a condition for the nontrivial solutions 
of the integral equation (fT4"]l to exist. This condition, 
which determines the frequency u, will be referred to as 
a dispersion equation. If the electron velocity is directed 
along the OX axis, the dispersion equation is 

[(1 - B S x )(l - B S Z ) - B 2 S 2 ] [1 - B S y ] = 0. (17) 

The dispersion equation (JTTj) describes all possible collec- 
tive oscillations of electrons and the dipole. If electrons 
are in the equilibrium state (vo = and Sq = 0), the 
dispersion equation can be decomposed into three equa- 
tions: (1-B S X ) = 0, (l-B Sy) = 0, and (1-B S Z ) = 0. 
The solutions of each of them correspond to different ori- 
entations of the induced dipole. It is evident that the so- 
lutions of equations (the frequencies of collective dipole 
and electron oscillations), which correspond to the x- 
and y-orientations of the dipole, are identical, because 
these two directions are physically equivalent in the ab- 
sence of charge carrier drift. The frequency of collective 
oscillations for the z-orientation of the dipole is, gener- 
ally speaking, different from those corresponding to the 
x- and y-orientations. 

In the general case, Eq. ([17]) is factorized into two 
separate equations. One of them is the equation 

(I - B Q S x )(l - B Q S Z ) = B 2 S 2 . 

One can easily verify that, for the frequencies, which are 
the solutions of this equation, the electric field of elec- 
trons has a symmetry, at which the induced dipole of a 
nanoparticle lies in the x — z plane, i.e. it corresponds 
to a mixed x ~ z orientation. The roots of this equation 
are two frequency branches. The branch, which is asso- 
ciated at vq — with the x-orientation of the dipole (i.e. 
which is determined by the equation BqS x = 1), will be 
called the x-branch. Accordingly, the other branch will 
be called the z-branch. As it was in the equilibrium case 
(vq = 0), the solutions of the equation 

BqS v = 1 
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correspond to the orientation of the dipole in parallel to 
the OY axis for any wo-value; therefore, the correspond- 
ing frequency branch will be referred to as the y-branch. 

The solutions of the dispersion equation (fT7| can be 
presented in the form 



A 



Rj(Vo,Qj,T p ), 



(18) 



where the subscript j ~ x,y, z denotes the x- 1 y- 1 and 
z-branches, respectively. We also used the notations 



Wpi 



Ulpl = 



2ne 2 no 



Vn = 



ojpih ' 



"0 — ; 1 — ; 



A = 



e 2 e* n o |(l|x|0}| 5 



The parameter A is responsible for the coupling between 
the 2DEG and the dipole; and Rj(Vo,flj,T p ) is a cer- 
tain complex-valued function depending on the electron 
drift velocity, the frequency of dipole oscillations, and 
the plasmon damping. For the x- and z-branches, the 
expressions for the functions R XiZ are 



Rx,z - 2 




(19) 



where the plus sign corresponds to the ir-branch, and 
the minus sign to the z-branch. It is easy to obtain a 
simpler expression for the y-branch, 

Ry = S y ■ 

The real part of this function, R'j , describes a shift of the 
oscillation frequency of the system with respect to the 
eigenfrequency of the dipole owing to the interaction be- 
tween the dipole and electrons. The corresponding imag- 
inary part, R" , is responsible for the additional damping 
(or growing) of oscillations in the hybrid system. 

From the analysis of integrals (f!6| . it follows that, 
in the extremely high-frequency case (j4]), the functions 
Rj practically do not depend on the parameter T p . (In 
the hydrodynamic approximation, which is used in this 
work, the plasmon damping is exclusively governed by 
the one-particle relaxation of electrons. The Landau 
damping for plasmons, which arises, if the electron kinet- 
ics is considered 22j . is assumed small.) In the absence 



of electron drift, it can be shown analytically that, if 



0, 



Rx v = VP I dk 



R = 2R X 



h,3 p 2/c 



(20) 



The nonzero value of i?" and, respectively, the oscillation 
damping (analogous to the Landau damping (22j) in the 
hybrid system originate from the collective interaction. 
Really, since the spectrum of plasmons (Eq. (fTU)) ') is con- 
tinuous, for any given frequency luq, there can always be 
found "resonance" plasmons, the charge waves of which 
are proportional to exp[i(k r r — tut)], where the wave vec- 
tor k r satisfies the condition w ± (k r ) — u>o- At the same 
time, the electric field of the dipole is a sum of Fourier 
components Ed.k<3xp[ikr — iujQt], with k = k r inclusive. 
The dipole field waves with the wave vectors k r and the 
"resonance" plasmons are cophased, and they propagate 
with the same velocity. This means that the correspond- 
ing dipole field acts permanently (without variation in 
time) on "resonance" plasmon charges. As a result, the 
work of the dipole field over the charges is nonzero, and 
the dipole loses its energy. Under nonequilibrium condi- 
tions, the dipole can acquire the energy from the electron 
subsystem. A necessary condition for that is 



R^(V ,Q 3 ;T P ) >0. 



(21) 



The sufficient condition for the hybrid system to be un- 
stable in whole is more rigorous: 



A_ 



R , '(V ,il j ,T p )>T , 



(22) 



If the inverse inequality is obeyed, either the system is 
stable or the oscillations in the system attenuate. 

4. Interaction between a Shallow Donor and 
Drifting Two-Dimensional Electrons 

In this section, we apply the results obtained above to a 
specific hybrid system consisting of a shallow impurity 
center and a heterostructure with a quantum well. 

It is well known that one-particle Coulomb impurities 
in semiconductors are characterized by low binding en- 
ergies, and the allowed photo-induced dipole transitions 
between impurity states correspond to the THz spectral 
range. Such impurities can be regarded as hydrogen- 
like atoms, the energy spectrum and the wave func- 
tions of which can be calculated in the effective mass 
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Fig. 2. Dependences of the real (a) and imaginary (ft) parts of the function Rj on the electron drift velocity at a fixed frequency of 
dipole oscillations (Qq ra 0.58). The other parameters are given in the text 



approximation (see, e.g., work [23|). The energy differ- 
ence between the basic, IS*, and excited, 2P, states is 
evaluated as E$-p — 3e 2 /8Kas, where as = kH 2 /me 2 
is the radius in the ground state. For example, for GaAs 
with to = 0.067mo and k — 12.9 (too is the free electron 
mass), we obtain E$-p ~ 4.12 meV and ap ~ 10 nm. 
This energy difference corresponds to the frequency 
uj w 6.2 x 10 12 s" 1 (» 0.99 THz). The phototransition 
S f-> P is an allowed electric dipole transition with the 
transition matrix element (l|x|0) « 0.52a|j. The cited 
parameters allow the polarizability of a Coulomb impu- 
rity to be calculated using relation (J7|)- 



To achieve conditions needed for the instability exci- 
tation, a substance for the heterostructure must be se- 
lected, which would be characterized by high electron 
velocities. Consider a quantum well on the basis of 
InAs with GaAs-barriers. It is known [24] that the effec- 
tive electron mass in InAs is small, m ps 0.023too, and 
the electron mobility is high even at room temperature, 
/i«8x 10 4 cm 2 /(V x s). It gives rise to very high drift 
velocities of electrons, up to vq » 6 x 10 7 cm/s [2J,|25|. 
In quantum wells on the basis of InGaAs, drift veloci- 
ties of the same order are observed. The difference be- 
tween the dielectric constants of the quantum well and 
the barrier can be neglected [20]. For numerical calcula- 
tions, let us choose such physical parameters that crite- 
ria (QJ are satisfied. In particular, let us fix the electron 



concentration no = 10 cm and the distance from the 
2DEG to the donor h = 4 x 10 -6 cm. The corresponding 
characteristic parameters, which were introduced by re- 
lations (HU), are tu p i 1.07 x 10 13 s" 1 , A « 0.0013, and 



r p w 0.03. The drift velocity of charge carriers is nor- 
malized by the quantity u p ih « 4.28 x 10 7 cm/s. 

In Figs. 2 and 3, the dependences of the real, Rj, 
and imaginary, R'L parts of the function Rj on the 
normalized velocity Vq and frequency f2 , respectively, 
are depicted. Figures 2, a and 3, a bring us to a con- 
clusion that, for the selected fixed parameters, the fre- 
quency of oscillations in the system is a little (in com- 
parison with wo) shifted toward lower values in the fre- 
quency x- and y-branches, and to both sides in the z- 
branch. Figure 2 also demonstrates that i?" changes its 
sign at Vb ~ 0.8, which corresponds to the drift velocity 
vo ps 3.4 x 10 7 cm/s. If condition (|2"2")l is also obeyed at 
that, an instability must emerge, and oscillations in the 
system start to grow. The figure also illustrates that the 
x- and y-branches merge together at low enough drift ve- 
locities, which is associated with the equivalence of those 
directions at Vq — > 0. Figure 3 testifies that the lower the 
frequency of dipole oscillations, the larger is the insta- 
bility increment. As to the y- and z-branches, R' y ' < 
and R" < in the considered ranges of frequencies and 
drift velocities; therefore, the corresponding oscillations 
in the system attenuate here. 

As was shown in works the inverse lifetime 

of donor electrons at shallow donors in GaAs can be of 
the order of 10 7 s _1 . For V « 0.94 and il Q « 0.58, we 
obtain --^R"ujpi « 5.28 x 10 9 s^ 1 . In this case, criterion 



is satisfied well, and the indicated instability effects 
must be observed for the hybrid system concerned. Note 
that the increase in the electron concentration or the 
distance h leads to a reduction of the instability incre- 



ment, because A 



_1 and A cx h 2 . We also verified 
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Fig. 3. Dependences of the real (a) and imaginary (b) parts of the function Rj on the frequency f2o at the fixed electron drift velocity 
Vo ~ 0.94. The other parameters are given in the text 



that all effects discussed depend weakly on the parame- 
ter T p . 

Therefore, using an InAs quantum well and a shallow 
hydrogen-like donor in the GaAs barrier as an exam- 
ple, we showed that an instability can take place in the 
a;-branch of collective oscillations in the system. Oscil- 
lations corresponding to the y and z frequency branches 
attenuate. 

5. Charge Waves 

In the previous sections, the dispersion equation for col- 
lective oscillations of the dipole and the 2DEG was de- 
rived and analyzed. Its solutions are the characteristic 
frequencies (eigenvalues) of oscillations. In this section, 
we are going to construct eigenfunctions, which corre- 
spond to those eigenvalues and are the solutions of the 
integral equation (fT4")h The analysis of these solutions 
will allow one to understand the behavior of the elec- 
tron subsystem at collective oscillations. 

Let us recall that the integral equation (| 14[) for the 
charge carrier concentration rif. must be specified for the 
corresponding type of collective oscillations, after the so- 
lutions of the system of algebraic equations (fT5|) for I x , 
I y , and I z (they are functionals of rife) have been found. 
For two mixed x — z orientations of the dipole, the quan- 
tity I y is nullified. Then, using the first and third equa- 
tions of this system, it is easy to obtain 



I, 



KZ X I X 



where the notations K x = BoSo/(l — BqS z ) and 
Kz = BqSq/(1 — BqS x ) are introduced (in these nota- 
tions, the dispersion equation for the x- and z-branches 



reads K X K Z = 1). Then, the solution of the linear inte- 
gral equation (|14[) . which corresponds to the x-branch 
of frequencies, is 



h ^ °-°o — 7—7 n — e 

A e (uj x ,k) 



(23) 



where the superscript in parentheses in the notation 
means that the frequency entering into the quantities 
Bo, K x , and A e belongs to the cc-branch. Generally 
speaking, the constant C is an arbitrary complex value 
characterized by an amplitude and a phase. The latter is 
insignificant, because it can always be zeroed by shifting 
the start of time counting. Therefore, C is assumed real- 
valued below. 

Substituting Eq. into the formula of the Fourier 
transformation, we obtain the space-time distribution of 
a concentration perturbation, 



n {x) (r,t) = CB Q J d 2 k 



k (k x + K x k) e -kh+ikr-iu x t 

A e (uj x ,k) 



(24) 



It is the real part of expression (|24[) that has a physical 
meaning. Let 



J x (r) = CB Q f d 2 k- 

J A e (w,-,fc) 

J z (r) = CB f d 2 k f t 
J A e (wj,fc) 



-A,'/i+ikr 



and Tj = u>jt, where j — x for the x-branch. Extracting 
the real part of the concentration perturbation (|24[) . we 
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15. 

2 



Fig. 4. Contour plots of the space-time dependences of electron 
concentration perturbations at Vb = 0, which correspond to the 
x-branch of frequencies: (a) the spatial dependence at r x 
(ft) the space-time dependence at Y = 

obtain 



(J' x (v) + K'Mt) - K»J f J(r)) co S r x - 



+(4'(r) + K' x J'J(r) + K' x 'J' z (r)) smr x 



(25) 



Solutions (|2"51 have the following properties. The func- 
tion n'( x >(r, t x ), as well as n^ x '(r,t), is even with re- 
spect to y. Since — f- <C 1, the amplitude of this func- 
tion changes slightly within several periods. At Vb = 0, 
n'( x \r, t x ) is an odd function of the coordinate x. In gen- 
eral, n'( x >(r, t x ) describes the behavior of charge waves 
in the two-dimensional space and in time. 

In Fig. 4, the contour plots of the space-time de- 
pendences of electron concentration perturbations cor- 
responding to the x-branch (in this case, the dipole is 




-2 2 



Fig. 5. The same as in Fig. 4, but at Vb ~ 0.94 

oriented along the OX axis) in the absence of electron 
drift (Vb — 0) are exhibited. In particular, in Fig. 4, a, 
the spatial distribution of the concentration at the time 
moment t x = 3ir/2 is shown. It has a minimum and a 
maximum to the left and to the right, respectively, from 
the coordinate origin, where the dipole is located. This 
distribution arose owing to the corresponding arrange- 
ment of dipole charges at that moment. The time-space 
illustration (Fig. 4,6) demonstrates that those maxima 
and minima alternate in time with a period of 7r. The 
spatial scale of perturbations has an order of several h. 

The contour plots of the space-time dependences of 
the electron concentration in the case where electrons 
drift along the OX axis (Vb ~ 0.94), are shown in Fig. 5. 
In the general case, if electrons drift, the dipole be- 
comes oriented in the plane y = (see the next section) . 
One can see that the concentration minima and maxima, 
which are exhibited in Fig. 4, are "blown" by the electron 
stream, with the perturbation boundaries moving faster 
than its center. The nonzero phase velocity of plasmon 
waves, which is perpendicular to the electron drift, re- 
sults in a smearing of perturbations located far from the 
dipole. 

A similar analysis of the time-space distributions of 
electrons can be carried out for the z-branch of the fre- 
quency dispersion. For this case, instead of Eq. (f2~5j) . we 
obtain 



(r. 
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Fig. 6. The same as in Fig. 4, but for the z-branch 



(K' z J' x (r) - K'JJ'J(t) + J' z (r)) cosr 2 + 




Fig. 7. The same as in Fig. 5, but for the z-branch 

somewhat more complicated than that obtained above 
for the x-branch. Perturbations located near the dipole 
move against the electron stream, whereas the remote 
perturbations move along the stream. At Jw 3.5, the 
character of perturbations changes, and the role of the 
wave propagation becomes dominant. 

Consider the solution of the integral equation (f!4|) for 
the dipole y-orientation. The real part of an electron 
concentration perturbation looks like 



-(K' z J':(r)+K':j' z (r) + J'J(r)) S inr z 



Ay) 



(26) 



The properties of this function are a little similar to those 
of the function n'^^r, t x ). In particular, this function 
is also even with respect to y. However, in contrast to 
n'( x \r, t x ), n'^'(r, r z ) is an even function of the coordi- 
nate x at Vo = 0. 

Figure 6 illustrates the behavior of the electron con- 
centration for the z-branch of frequencies at Vq — (the 
dipole is oriented along the OZ axis). The spatial dis- 
tribution of the electron concentration calculated for the 
time moment t x = 3tt/2 (see Fig. 6, a) has a maximum, 
which is located under the dipole and is symmetric with 
respect to the substitutions x — > —x and y — > —y. Fig- 
ure 6,& demonstrates the behavior of a perturbation in 
time, which is close to periodic. Figure 7 exhibits the 
same as in Fig. 6, but provided that electrons drift with 
the velocity Vq w 0.94. In this case, the dipole is also 
oriented in the plane y = 0, but the wave dynamics is 



{v,T v ) = e-'y Tv [j£ (r) cos r„ + Jesuit,,], (27) 



where the notation 
J y (r) = CB Q J d 2 k 



A e (u> y , k) 



V g-fc/i+ikr 



is used. The properties of n'^(r,T y ) wee also simi- 
lar with to those of n'^ (r, t x ). At Vq — 0, the fre- 
quency branches coincide; therefore, in order to de- 
rive n'( y )(r, T y ) from n'^(r, t x ), the substitution X O Y 
should be made. The same concerns Fig. 4 as well. 
In contrast to n'^(r, t x ), n'^ v '(r,T y ) is an odd func- 
tion of the coordinate y for an arbitrary Vo, which 
is connected with a constant orientation of the dipole 
in parallel to the OY axis. Figure 8 illustrates the 
coordinate (panel a) and coordinate-time (panel b) 
dependences of the electron concentration perturba- 
tion, which corresponds to the y-branch at the elec- 
tron drift. At Y = 0, the integrand in the integral 
Jj^r) is an odd function, so that n'^ v \x,y — 0,r y ) = 0. 
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Fig. 8. The same as in Fig. 5, but for the y-branch at Y = 0.5 



Therefore, Fig. 8 was plotted for y = 0.5. It can 
be interpreted like Fig. 7,6: perturbations are blown 
down the electron streams; however, in this case, per- 
turbation minima transform into maxima and vice 
versa. 

Hence, in this section, we determined the behavior 
of perturbations in the two-dimensional electron gas 
in the hybrid system for various types of solutions, 
which takes into account whether or not the electrons 
drift. 



6. Behavior of the Dipole Induced in a 
Nanoparticle 

Consider the behavior of the nanoparticle polarization 
during collective electric oscillations in the hybrid sys- 
tem. This behavior is described by the dynamics of the 
induced dipole moment ([5]) . To determine the dipole mo- 
ment, we have to calculate the electric field created by 
electrons at the dipole localization point making use of 
definition ([9]), the first formula in system (fl"2)) . and the 
Fourier components of the concentration rik (the latter 
were analyzed for various dispersion branches in the pre- 
vious section). 

First, let us consider the solutions corresponding to 
the dispersion x-branch, provided that the components 



of the dipole moment are 

7 / \ . 27re/3 (w K ) x 
dx(Ux) = -tC = C x 

K 

d y (uj x ) = 0, 
4K) = iC x K x , 



(28) 



where C x is a new arbitrary constant. Whence, we ob- 
tain the law for the time evolution of dipole moment 
components, 



D x = e^ Tm (C' x cosT x + C x mnr x ), 
D z = -e^ T * \(C' X KZ + ClK' x ) cost, 
+ (Cx K " ~C' X K' X ) sinr x 



(29) 



System (1291) is a parametric equation of an ellipse: the 
end of the dipole moment vector moves along an elliptic 
trajectory in the plane y = 0. It is easy to verify that 
K' x > and K" > for the x-branch of the frequency 
dispersion; hence, this circulation is counter-clockwise. 
Since the dipole is in a self-consistent field with drifting 
electrons, the parameters of this ellipse depend on the 
drift velocity and the frequency of dipole oscillations. 

Let us put the arbitrary constant C x = 1 and, for sim- 
plicity, let exp (^7-t x ^J rs 1, i.e. we consider the behavior 
of the system within several periods. Then, 



(30) 



D x = cost,, 

D z = K' x sin t x — K'l cos t x . 



After excluding the parameter t x from system (|30|) . we 
obtain the equation for this ellipse in the form 



aDl + 2bD x D z + cDt = 1, 



(31) 



where a = 1+ {K^/K' x f, b = K'J/K' X \ and c = 

The rotation of the coordinate system by an angle a a 

determined by the equation 



1 — K' - K"~ 
tan Qa: = f—^- ± 



' 1 - K ' - K" 

X X 

2K" 



brings about the canonical ellipse equation, 



D 2 D 2 
A 2 B 2 



= 1, 



+ 1, 
(32) 

(33) 



where D Xg and D Zo are the dipole components in the new 
coordinate system, and A and B are the principal ellipse 
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axes. 



A 2 = 



1 + tg 2 a x 



2btga x + c tg 2 a a 



B 2 



1 + tg 2 a x 



atg 2 a x - 2btga x + c 



(34) 



The parameters A and B are used to introduce the ellipse 
eccentricity, 



B 2 , 



A 2 

1-^ *<B. 



(35) 



The e- magnitude varies from zero to 1. The ellipse 
transforms into a circle at e — > and into a line at e — > 1. 
The eccentricity of the ellipse, along which the dipole 
circulates, evidently depends on the electron drift veloc- 
ity and the frequency of dipole oscillations, because the 
principal axes depend on K x . 

Note that the time dependence of the dipole is very 
simple in the coordinate system coupled with the prin- 
cipal axes, namely, 



D Xo = Acost x , 
D Zo = Bsm.T x . 



(36) 



Similar formulas can also be derived for solutions as- 
sociated with the dispersion z-branch. In this case, 
the dipole also circulates along an ellipse, but now 
"clockwise". The principal axes of this ellipse depend 
on K z (ui z ). To obtain the canonical ellipse equation 
for the z-branch, the substitutions a <-> c, b —> —b, and 
K x {uj x ) — > —K z (u) z ) should be made in the equations 
written above. Now, the eccentricity is described by 
expression (|35[) . but the principal axes swap (A o B), 
and the angle a x changes to a z . The eccentric- 
ity is the same for both dispersion branches, because 

K x {n x ) ~ -k z (q z ). 

Concerning the y-orientation of the dipole, the eccen- 
tricity of the corresponding ellipse is equal to 1 for all Vo, 
which is associated with the dipole orientation strictly 
along the OY axis for any Vo- 

The dependences of eccentricity e on the electron drift 
velocity for three different dipole oscillation frequencies 
are depicted in Fig. 9, a. The solid curve exposes the ec- 
centricity that corresponds to the same frequency as in 
all previous figures, where this frequency was fixed, i.e. 
Qq ~ 0.58. Two other curves correspond to half as high 



and twice as high dipole oscillation frequencies. The fig- 
ure testifies that the eccentricity is equal to 1, if the 
drift velocity equals zero, i.e. the dipole oscillates with 
the frequency fi x along the OX axis and with the fre- 
quency Vt z along the OZ axis. Figure 9,6 shows that the 
lines, along which the dipole oscillates, are perpendicular 
to each other at Vq = 0. If Vo increases, the eccentricity 
decreases, and the line transforms into an ellipse. The el- 
lipse, which corresponds to the cc-branch, rotates "clock- 
wise"; therefore, the absolute value of the angle, by which 
its axes are rotated with respect to the OX axis, grows. 
The other ellipse rotates "counter-clockwise"; therefore, 
the absolute value of its orientation angle decreases. At 
a certain drift velocity, those ellipses coincide. 

A similar situation takes place at other dipole oscilla- 
tion frequencies as well. The curves in Fig. 9 also demon- 
strate that the eccentricities have a similar structure at 
different frequencies, and the smallest e- value is realized 
for a curve that corresponds to Hq » 0.58. At this fre- 
quency, the function Rj, the plot of which is exhibited 
in Fig. 2, changes drastically in a vicinity of the drift 
velocity, which corresponds to the eccentricity minimum 
(Vo ~ 0.43). The plots of the functions a x and a z inter- 
sect each other also in a vicinity of Vo w 0.43. At other 
frequencies, such eccentricity minima take place at dif- 
ferent drift velocities, and the behavior of the function 
Rj in vicinities of those minima (to be exact, in the re- 
gion, where the eccentricity diminishes sharply) is also 
drastic. 

Let us elucidate the physical reason for the nontriv- 
ial behavior demonstrated by the induced nanoparti- 
cle dipole. For this purpose, consider the electric field 
of plasmons at a certain point (x = y = 0, z = h) re- 
mote from the 2DEG. The dispersion law ([TBI has two 
branches. However, for simplicity, let us consider any 
of them, omitting the subscript at the frequency and 
neglecting the plasmon damping. Let a plasmon with 
the frequency uj(k) and the two-dimensional wave vector 
k = (k, 0) be excited. Then, at the distance h from the 
coordinate origin, it creates a field 



= — ik\Ao\e 
lole 



— kh—iujt 



E, 



E r = k\Ar ) i "- kh "-' lUt 



where \Aq \ is an arbitrary amplitude. Only the real part 
of the found field has a physical meaning. Its compo- 
nents are 

{E' x = -fc|A,|e- fc,l sin(wi), 
\E' Z = k\A a \e- kh cos(ujt). 

We see that, at a certain distance from the coordinate 
origin, the field created by plasmons circulates along a 
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Fig. 9. Dependences of (a) the ellipse eccentricity tj and (b) the rotation angle otj (in degree units) of the principal axes of the ellipse 
on the electron drift velocity 



circle. It is evident that the dipole moment of a nanopar- 
ticle, which is given by formula ([5]), when responding 
to an external force, will also circulate under the in- 
fluence of the plasmon field. When collective oscilla- 
tions take place in the hybrid system, not a single but 
many plasmons with different k are excited. The total 
field has a more complicated time dependence, which re- 
sults in the dipole circulation along an ellipse. Hence, 
the dynamics of the plasmon field is a key for under- 
standing the behavior of the nanoparticle polarizabil- 
ity. 

It is important that the behavior of the induced 
dipole in time can be experimentally observed owing 
to features of the emission by the hybrid system. Re- 
ally, having solved the problem of an electrostatically 
coupled nanoparticle and the two-dimensional electron 
gas, we determined the motion of charges in the hy- 
brid system. The motions of polarization charges in 
the nanoparticle and electrons in the quantum well 
correspond to an electric current depending on the 
coordinates and the time. The alternating current 
generates radiation. Let the current density be des- 
ignated as j(x,y,z,t). Then, the vector potential 
of the radiation field can be determined by the for- 
mula [53| 

ki js 1 f , j , i( x i,vuzi,t- f ) 

A{x,y,z,t) = — / dxxdyxdzi , (37) 

co J R 

where c is the speed of light, and R = [(xi — x) 2 + 
(yi - y f + ( z i - z) 2 ] 1 / 2 . Using formula J37} for the far 
radiation zone, we obtain the Fourier components of the 



vector potential in the form 

A(x,y,z,u j ) = 1 J> , (38) 

Hq 

where Ro = [x 2 + y 2 + (z — h) 2 ] 1 / 2 is the distance 
to the observation point, and fco = Wj/cq. The fre- 
quency ujj belongs to one of three frequency branches. 
The dipole moment, which determines A, is formed 
by the self-consistent system composed of the elec- 
trons and the nanoparticle. Hence, the radiation 
characteristics substantially depend on such parame- 
ters of the system as the distance between the dipole 
and the quantum well, the concentration of electrons, 
the frequency of dipole oscillations, and the electron 
drift velocity. At finite drift velocities of electrons, 
the induced dipole circulates along an ellipse, and 
this fact must be reflected in the radiation polariza- 
tion. 

7. Conclusions 

A hybrid system consisting of an isotropic nanoparticle 
and a heterostructure with a quantum well has been con- 
sidered. The nanoparticle is assumed to be polarizable 
in an external electric field and to have characteristic res- 
onance frequencies in the terahertz range. In this case, 
there emerge collective oscillations of the nanoparticle 
polarization and plasmons in the two-dimensional elec- 
tron gas of the hybrid system. 

Dispersion relations for the collective oscillation fre- 
quencies were obtained and analyzed. Possible frequency 
branches were determined and classified. Additional 
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damping of predicted oscillations was revealed, the origin 
of which is similar to that of Landau collective damping 
in plasma. The electron drift results in a reduction of the 
additional damping. At sufficiently high drift velocities, 
owing to the energy of the electric current, there emerges 
an instability in the system, and the oscillations in one 
of the dispersion branches grow in time. The electric 
instability increment increases, when the distance be- 
tween the dipole and the electrons diminishes and when 
the drift velocity increases. 

The predicted effects were illustrated using the re- 
sults of numerical calculations for a shallow hydrogen- 
like donor in the barrier of an InAs-based heterostructure 
with GaAs barriers, which were taken as an example. 
The space-time dependences of concentration perturba- 
tions in the two-dimensional electron gas in the course 
of collective oscillations were analyzed. The calculated 
behavior was demonstrated to be substantially different 
for different frequency branches, both in the absence and 
the presence of drifting electrons. 

The polarization oscillations of a nanoparticle were 
studied. It was found that the induced dipole is char- 
acterized by a complicated dynamics at nonzero drift 
velocities. In particular, in two of three branches, the 
dipole circulates along elliptic trajectories depending on 
the electron drift parameter. It was shown that the fea- 
tures in the nanoparticle polarization behavior could be 
observed by measuring the emission of the hybrid sys- 
tem. 

The practical interest to the new phenomena in hy- 
brid systems may consist in a capability to excite the 
emission by nanoparticles by applying an electric current 
and the electrically stimulated generation of THz radi- 
ation. These phenomena can also be used for the field- 
controlled addressing to individual nanoparticles, which 
is a key problem at the implementation of quantum cal- 
culations |28j . 
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y KBAHTOBIH 51MI 

B.O. KoueAan, CM. Kyxmapyn 

P e 3 k> m e 

Po3rji5myTO ri6pnflHi CHCTeMii, mp CKjiaflaiOTBca 3 HaHonacTimKH 
Ta HaniBnpoBiflHHKOBo'i reTepocTpyKTypn 3 KBaHTOBOio hmokx Ha- 
HonacTHHKa e TaKoio, mp nojispn3yeTBCH y CTopoHHBOMy ejieKTpn- 
HHOMy nojri. 06rpyHTOBaHO Ta ccpopiviyjiBOBaHO MOflejiB ri6pHflHoi' 

CHCTeMH. OTpiIMaHO TO"ffli p03B'5I3KII piBHHHB. 3HaHflem HaCTOTH 

KOjiiiBaHB 3apHfliB ri6pHflHoi' chctcmh Ta i'x ^o^aTKOBe 3aracaHHsi, 
mp 3yMOBjieHO B3aeMOfliero flnnojia 3 njia3MOHaivin. ilpnpofla flOfla- 
tkoboto 3aracaHHH nofli6Ha pp 3aracaHHH JIaHflay. IlpoaHajii30Ba- 
ho noBefliHKy b naci Ta npocTopi 36ypeHb KOHu,eHTpau,ii' flBOBHMip- 
hhx ejreKTpoHiB. ^focjiifljKeHO nojiHpH3amHHi KOjiiiBaHHs HaHona- 
cthhkh, SHaiifleHO, inp npn HeHyjiBOBHx flpencpOBHx hibh^koctstx 
HaBe^eHa nojiHpii3an,i5i xapaKTepH3yeTBC5i CKjiamroio miHaMiKOio. 
3oKpeMa, fljia flBOx i3 Tptox rijiOK nacTOTHoi' flHcnepci'f BeKTop no- 
jiapH3an;ii o6epTaeTBCH no ejiinranHnx TpaeKTopiax. y BnnaflKy, 
kojih flo KBaHTOBoi hmh npHKjiafleHe nojie Ta Tene CTpyin, 3ara- 

CaHHH 3MiHK36TBCH Ha HapOCTaHHH KOJIHBaHb ri6pHflH01 CHCTeMH 

y naci, inp BimTOBi^ae ejieKTpHMHiH HecTiHKOCTi ri6pHflHoi' CHCTe- 

MH. HOBI HBHina B ri6pH^HHX CHCTeMax MO^KyTB 6yTH 3aCTOCOBaHl 
flJIH 36yflJKeHHH BHnpOMlHIOBaHHH HaHOHaCTHHOK CTpyMOM Ta fljia 

ejieKTpHHHoi reHepani'i BHnpoMimoBaHHH b TeparepnoBift o6jiacTi 
cneKTpa. 



B3A6MO/1I51 I30TPOI1HOI HAHOHACTHHKH 
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